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Abstract 

We use semi-classical and perturbation methods to establish the quan- 
tum theory of the Neumann model, and explain the features observed in 
previous numerical computations. 
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1 Introduction 



In a previous publication pQ, we have explained how to determine numerically 
the spectrum of the quantum Neumann model. The Neumann model describes 
a particle on a sphere subject to a quadratic potential. It has been introduced 
in the nineteenth century by Neumann [2J, who showed that it is integrable. 
More recently, Moser [3] and Mumford have shown how it fits in the modern 
approaches to algebraically integrable systems and generalized the model to any 
number of dimensions. More to our point, the quantum version of the model 
is also integrable jSE]. This means that its solution reduces to the resolution 
of a one dimensional Schrodinger equation, as it is the generic behaviour for 
algebraically integrable systems. The notable feature compared to the textbook 
examples is that the different conserved quantities, which all appear in the same 
equation, must be solved for simultaneously. This is a non trivial problem which 
has been attacked in p]. We will precise and prove some statements made in this 
paper. Moreover, the Neumann model is related with Stackel systems as noted 
by Gurarie and fully developed recently in jSj. 

We will start with a brief presentation of the model and the equations which 
have to be solved. The physics of the model is driven by a general scale pa- 
rameter v, which is proportional both to the strength of the interaction and the 
radius of the sphere. Small values of v correspond to a quantum regime and 
are next studied by perturbation theory. Then, we introduce the semi-classical 
analysis which is appropriate for large values of v. The comparison of energy 
levels obtained by semi-classical methods and by the numerical solution of the 
Schrodinger equation shows that it is indeed the case. However, we find a small 
constant discrepancy between the two. In a last part, we recover this constant 
by analytic studies of the model at large v according to these two methods. 



The Hamiltonian of the system is given in terms of angular momenta J^i = 



with the constraint ^x\ = r 2 . We assume a\ < a 2 < • ■ ■ < ajv- One can always 
shift all the a>k by a constant since this merely adds a constant to H. The system 
is classically and quantum mechanically integrable since the following quantities 
commute, classically and quantically: 



2 The model. 
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To solve the equations of motions, one introduces the "separated" variables 
tk given as the roots of: 



Positivity of x\ implies that there is exactly one root in each interval [a*., 
It has been shown in |H| that finding the common eigenvector for the F k with 
eigenvalue /& reduces to solving the unique one-dimensional Schrodinger equa- 
tion: 



This is a linear differential equation with singularities at the a k and at infinity. 
The situation is explained in more details in £Q. 

It can be useful to write the potential term in the form: 



with v = J2k fk = "l2k x k = r<1 ■ Classically the action S satisfies the separated 
Hamilton- Jacobi equation: 



In Figure 1 we draw the different possibilities for the potential in the Neumann 
case N = 3, with oscillator strengths 0, 1 and y. 

For v large, we are in case A), b\ is close to 1 and b 2 is close to y. Classically 
the allowed region is such that for each t k one has P(t) = Yl(t~ a n) Y\ii — b m ) < 
so that, since t k e [a k , a k +i\ one has b\ < t\ < 1 and b 2 < t 2 < y. As v is lowered, 
both b\ and 6 2 move to the left, and when v — > 0, one has 61 — » —00 so that one 
is in case C) or D). One goes from case D) to case A) either through the cases B) 
or C). In case B), b\ and b 2 both belong to the interval [0, 1] so that the allowed 
region corresponds to bx<t\< b 2 and t 2 describes the whole interval [l,y]- In 
case C) and D), b\ is negative and b 2 is on either side of 1. One of the t k describes 
the whole interval and the other is constrained away from 1 by b 2 . It is easy to 
check that no other configuration is compatible with the existence of t\ and t 2 
satisfying P(t) < in the appropriate intervals. In the general case, we will show 
in the sequel that all the 6j are real and belong to the interval [a$_i, Oi+i], with 
a supposed to be —00. There is therefore a total of 2 N ~ 1 possibilities for the 
relative positions of the 6j and a k . 

At v — 0, the system reduces to a free particle on the sphere. The energy 
levels are given by the usual formula + 1) for spherical harmonics, but 
with a different choice of basis in the (2j + l)-dimensional eigenspace. These 
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Figure 1: The four different possibilities for t 1? t 2? with y = 2. 



basic functions are called spheroidal harmonics. These harmonics diagonalize the 
perturbation by the potential, so that energy levels separate linearly for small v. 

For large v, the potential term becomes dominant and the point is restricted 
to the vicinity of the two opposite poles x% = ±r, x 2 = 0, x 3 = 0. In the 
next approximation, we can neglect the curvature of the sphere and we have two 
independent harmonic oscillators in X2 and x%. The configuration around the two 
poles are linked by an exponentially small tunnelling amplitude, so that all states 
come in almost degenerate pairs. 

3 Small v perturbation theory. 

In the case v = 0, the singularity at infinity of the equation (J3J) becomes reg- 
ular. Even if the behaviour of the solution at infinity has no direct bearing on 
the physical wave function, it appears important for the characterization of the 
solutions and provides the quantization condition for the energy. 

Indeed, recall that around a regular singularity at x = 0, a second order 
differential equation has two independent solutions of the form x a (l + bx + • • •), 
with the exponents a roots of a second order equation. At the finite singularities 
a*;, the exponents are and 1/2 with respective monodromies +1 and —1. In our 
previous work [T], we determined that to produce a single valued wave function 
on the sphere, the solutions must be of definite monodromy at each of the points 
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a*;. Monodromy at infinity is the product of the monodromies at the a& and is 
therefore ±1. Since the exponent at infinity is solution of: 

and a must be integer or half integer to get the right monodromy, we conclude 
that the energy E is ^h 2 j(j + N — 2) with any non negative integer j = 2a. 

For a given energy E, there is a second root a which yields a solution vanish- 
ing at infinity. However it cannot be the solution we are seeking with diagonal 
monodromy at the a*., since such a solution would be single valued and regular 
on a compact Riemann surface. In the case where all monodromies are one, the 
solution has just one singularity, a pole at infinity and is therefore a polynomial of 
degree a. These functions are known as spheroidal harmonics of first species Jill . 

In order to find these polynomial solutions, we write them in the form FJ(t — 
Xj). Inserting in equation Q, we get: 

E (t _ A . )(t _ A .) + 2^ (t - a k )(t - \j) "4/^r^ = ° (6) 

As it should be, this equation requires that fk = and we recover the value 
for the energy from the dominant term at t — > oo. This equation is equivalent to 
the condition that its residues at the ak and Xj all vanish. The residues at the Xj 
give a system of non-linear equations which determine the A.,: 



The residues at the ak then determine the conserved quantities as: 

3 J 

The system of equations (J7J) can be seen as the condition of equilibrium of a 
system of points with logarithmic repulsive potentials, some of which being fixed. 
In this formulation, it is therefore clear that there is a unique solution with a 
definite number of Xj in each of the intervals [ak, flfc+i]- In the case N = 3, this 
gives d + 1 different solutions of degree d. Similar calculations can be done for 
solutions which are of the form Yl kEl Vt ~ a k Ylj(t ~ ^i) wri ere J is a subset of 
the singularities. This gives spheroidal harmonics of the second, third and fourth 
species jH] for N = 3. The total number of these harmonics with energy + 1) 
is exactly 2j + 1, the known degeneracy of this energy level. 

Numerical solutions of the system of equations (JJJ) are easy to get by a mul- 
tidimensional Newton method and have been used to produce initial solutions in 
our previous paper. 
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From the knowledge of this v = solutions, we will describe the behaviour of 
the solutions for small v. As soon as v is different from zero, the solution gets an 
essential singularity at infinity. We therefore multiply the spheroidal harmonic 
by the simplest function singular at infinity, namely exp(pt). In the vicinity of 
v = 0, we expect p to be small and neglect higher powers of p. Inserting in 
equation (jSJ) gives: 

\ ^ 1 1 \ ^ 1 J_ \ ^ fk 

^(t- \ t )(t - Xj) 2 2^ (t - a k )(t - Xj) m^t- a k 



1 1 1 



j 



The term proportional to \jt in this equation immediately gives that p is pro- 
portional to v, so that the hypothesis that p is small is valid. 

2np + ^Np - -^v = 

Here n is the degree of the polynomial. The roots Xj of the wave function satisfy 
equations which are small perturbations of equation 0, hence are close to their 
values for v = 0. The term in 1/t 2 then gives the perturbation of the energy: 



E = h 2 n(2n + N - 2) + v 



2(N + 4n) 



Similar calculations can be done for the spheroidal harmonics of the other species 
and this nicely reproduces the slopes of the numerical results at small v. 



4 Semiclassical analysis. 

As an integrable system, the Neumann model lends itself to a simple semiclassi- 
cal analysis writing basically that the Bohr-Sommerfeld conditions quantize the 
Liouville tori. Equivalently, the separated Schrodinger equation can be analyzed 
by the WKB method. 

4.1 The WKB treatment. 

The presence of a singular factor in the first derivative term of the equation (JHJ) 
seems to be a problem for the application of the WKB method in our case, but we 
will show that in contrast to what happens in the spherical symmetry case 
this does not cause troubles. In a first approach, we write an expansion in powers 
of h for the action S such that ^ = exp(iS/h): 

s = s + hs 1 + h 2 s 2 + ■■■ (9) 
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As usual the dominant term in h 2 does not depend on the first derivative term 
in the equation and gives for S : 



S = dt 



4ll(*-afc) 
The following term gives: 

is* - 2S>s[ + hs[ t 1 - = o 

Due to the form of So, the singular terms at cancel and one gets: 

exp(iSi) = Y[(t - bj)~* 

At this order, the wave function is regular at the and one can also show that 
the following correction to S is regular at a^. In contrast, the first correction is 
mildly singular at the turning point b\ and the following correction S2 has so bad 
singularities that the expansion breaks down. In the cases where the b\ are away 
from the a^, the usual connection formulae around the turning point ^0] justify 
an additional 7r/4 phase factor in the quantization conditions. 

Another approach is to desingularize the equation (jSJ) by introducing a new 
variable x defined by the (hyper-)elliptic integral [§]: 

x= / - (10) 



VTI(*-a*) 
The application of the chain rule gives: 

d 2 m 1^ 1 dm _ 1 d 2 ^ 
lit 2 + 2^t-a k ~dl ~ U(t - a k ) dx 2 

k 

so that the Schrodinger equation becomes: 

d 2 m v i — r . , . . . 

^-i^nc-«*=o (id 

This equation has neither singularities nor first derivative term, so that the usual 
WKB formula gives: 

9{x)* (n^-^exp Q J dx^-Hit-b^ 

It should be emphasized that the variation of x is alternatively real and purely 
imaginary when t crosses the a k , introducing additional i factors in the action. 
In the traditional iV = 3 case, x is an elliptic integral and t can be expressed 
with Weierstrafi elliptic functions. When t goes from —00 to +00, x describes 
the boundary of the half-periods rectangle. In the cases where there are more 
singularities, it is possible that the t variable is not well defined as an inverse 
function of x, however this is of no consequence for the study of the Schrodinger 
equation which proceeds separately on each of the intervals [ak,dk+i]- 
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4.2 Remarks on the zeroes of the potential. 



The above elliptic formulation allows to demonstrate quickly the claim we have 
made that there must be in each interval [oj, a region where the product 
P(t) = Y[{t — a n ) n(^ — b m ) takes negative values. Recall from pQ that the wave 
function is either even or odd under each parity operation Xk — > — In terms 
of the t variable, this translates into monodromy around a^. Now, the x variable 
of the elliptic parametrization behaves as x(a,k) + c^y/t — % in the vicinity of a& 
so that the wave function is even or odd as a function of x — x(ak), which means 
that either \l/ or its derivative vanishes at each of the x{ak)- 

Multiplying eq. by \1/ and integrating by parts, we then get for each k: 



Under the elliptic parametrization, the image of the real axis is a sequence of 
alternatively real and purely imaginary intervals, according to the sign of the 
product Y[{t ~~ a k)- The ty' 2 (x) term has therefore the same sign as this product. 
Since \l/ is real, Eq. (|T2*j) cannot be satisfied if the potential remains of the same 
sign as the kinetic term. This is equivalent to the condition P(t) < at some 
point of the interval [a,, aj + i]. Finally, we prove that this allows to establish 
that quantum mechanically the bj obey the constraints we have discussed in the 
framework of the classical analysis. 

The relative position of the and bj depends on the signs of the The 
only one with a fixed sign is J'n, which is positive. All combinations of the signs 
of the other ones are possible, which gives 2 N ~ 1 possibilities. We will show that 
each of these sign combinations corresponds to a unique set of relative positions 
of the bj. The fact that we identify N — 1 necessary real zeroes of the potential 
shows that all bj must be real. 

More precisely, we will show that bj belongs to the interval ]a,j, dj+i[ if fj is 
positive and to the interval ] <2j_i, CLj \ otherwise. This means that we have to show 
that in the interval ]a,j, dj+i[, there is exactly one zero if fj and fj + \ are of the 
same sign and two if fj is positive and fj + \ is negative. If the / are of the same 
sign, the potential reaches opposite infinite limits at the ends of the interval, so 
that it must have a zero in between. In the other case, the potential goes to 
positive infinity at the two ends of the interval and we proved that the potential 
must be negative somewhere in the interval. This is only possible if the potential 
has two zeroes which we identify as bj and bj+\. The only remaining thing to 
prove is that there is a zero in the interval ] — oo, ai[ when fx is negative. In this 
case, we use the fact that fk is positive, so that the dominant term near — oo is 
negative. The potential is positive at 0— hence there must be a zero in between. 
Since we identified N — 1 zeroes, we have exhausted the whole set of zeroes for 
the potential and there cannot be any other, proving in particular that all zeroes 
are real. 




(12) 
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4.3 The quantization conditions. 

In fact, using any of the two WKB approaches, we get the same semiclassical 
quantization conditions: 

f as+1 , / v not - bi) i, „ 

in the case where all the monodromies are one and bj is in the interval [a,j, 
The integer rij is the number of zeroes of the wave function in the interval 
[a,j , a-j+i] ■ If the monodromy around dj+i is —1, we should add one half to rij. 
In the cases where bj is near a,j or becomes smaller, the 7r/4 phase factor is no 
more a good approximation and we will postpone to further studies the relevant 
phase factors. This equation has already appeared in the literature, for example 
in |E1IZ|, where the refinement with the ir/A phase first appears in the last paper. 
It is introduced there using the machinery of Maslov indices. In this case, the 
justification through the connection formulae is very direct, it simply uses the 
linear approximation of the potential around the turning point and consequently 
the asymptotic analysis of the Airy function: 

Ai{x) ^ ^jc-f^, M{-x) ~ +oc cos(^ 3 / 2 - \) 

In order to compare the semiclassical predictions to the exact numerical solutions 
found in fl., we solve numerically the quantization conditions (fl"3|) and plot the 
resulting expression for the energy as a function of v compared with two exact 
energy levels which are degenerate in the large v limit. We do this for a small 
number of energy levels in order to keep the graphs uncluttered and plot them for 
oscillator strengths equal to (0, 1, 2) in Fig. 2 and (0, 1, 5) in Fig. 2, both in the 
case N = 3. In either figure, the three heavy lines are the WKB predictions. We 
see that the small v regime is poorly accounted for, which is of no surprise since 
in this case, we are no longer in the situation A) of Fig. 1 and their should be 
phase corrections different from it/ '4 and moreover, we are in the deep quantum 
regime, where we do not expect WKB approximation to be good. Nevertheless, 
we plan to study in further works the different cases of Fig. 1 and the transitions 
between them. 

The most notable feature of the curves plotted in Fig. 2 and 3 is that in the 
large v limit, there is a small constant separation between the exact and WKB 
curves. This separation is independent both on the energy level and the oscillator 
strengths. This contradicts the naive expectation that the corrections to the 
WKB approximation are of order h 2 with respect to the order zero solution which 
has an energy of order 1/h. This order of the correction is true for the solution 
away from the turning points, but the difference from the Airy function may 
create order h correction to the phase as noted in JUj- It is however remarkable 
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Figure 2: Energy levels for y = 2. 

that this finite correction to the energy does not depend on the level, so that the 
transition energies have vanishing errors in the small h regime. It is the purpose 
of the following section to show through perturbative computations the validity 
of this result. 

5 Large v study. 

We first compute perturbatively for large v or equivalently small h the energy 
levels as obtained from the "enhanced" Bohr-Sommerfeld rules, and then as 
perturbations of quantum oscillators sitting around the pole of the sphere. 

5.1 Semi— classical calculation. 

We limit ourselves to the N = 3 case and we fix the oscillator strengths to (0, 1, y) 
as in pp. Since v is large, the point on the sphere is limited to the vicinity of 
the points (±1,0,0), whose separated variables are t\ = 1 and t 2 = y. With b\ 
close to 1 and b 2 close to y, t\ is classically limited to the interval [b\, 1] and t 2 to 
the interval [b 2 ,y}, which limits the point to the vicinity of the pole. To do our 
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Figure 3: Energy levels for y = 5. 

perturbative calculations, we will write b\ = 1 — a and hi — U ~ (3, so that a and 
/3 are small. The integrand in the Bohr-Sommerfeld rules takes the form: 



(t-l + a)(t-y + {3) 



t(t-l)(t-y) 

When integrating in [bx, 1], we cannot expand the factors t — l + a and t—1, but 
all the others can be expanded around t — 1 and in powers of j3. The resulting 
integrals are easily computable yielding: 

1 1 2 1 a/3 

- an H a Ti n + ■ ■ ■ 

2 16 4j/-l 

Similarly for the integral in [62, 2/] we keep the factors t — y + [3 and t — y and 
expand the others around t = y in powers of a. We then get the integral: 

1/3 1 (3 2 1 a (3 

7T H -— 7T H ; 7T + • • • 

2 16 2/3/2 4 ^ _ !) 
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We write that the first integral is equal to ^i={n\ + |) and similarly the second 
equal to ^=(n 2 + \) which determines a and (3: 

Note that this justifies a posteriori the claim we made that a and /3 are small, and 
become smaller when v increases. Finally the energy E = \v{l + y — b\ — 62) — 
|v(a + f3) takes the simple form, in which there is no mixing between n\ and n 2 
excitations: 

E WKB = Km + -) + ^{n 2 + -)J 27 27 (15) 

Of course the first term corresponds to the energy of two independent harmonic 
oscillators of respective excitations n\ and n 2 , and nicely fits the ^/v look of the 
energy levels at large v apparent in the figures of jT| . The second and third term 
shift these square root terms towards their correct position, but miss it by a small 
constant. As we show below this small constant is accessible to a more correct 
quantum computation. 



5.2 Quantum calculation. 

The principle of this calculation is to work out a perturbation of harmonic os- 
cillators. We detail it in the case N = 3 and it is then convenient to go to the 
elliptic parametrization given by eq. (|10p. We expand the change of variable in 
the vicinity of t = 1 and t — y. Taking the origin in t = y, we write 

/dr 1 / r r 
T^vW 3 ^ + 2^ + 2(y-l) + ' 

which can be inverted to yield 

r = y-t = l -y(y - l)x 2 - ^y(y - l)(2y - l)x A + ■■■ (16) 

A similar computation gives the parametrization around t = 1 

r = l-t=- A {y- l)x 2 - l(y - l)(y - 2)x 4 + ■ ■ • (17) 

These values of the variable t will be used in the Schrodinger equation ([11)1 to 
express in terms of the variable x the potential ±(v/4h 2 )(t — — b 2 ). The sign 
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can change since in the calculation of x, we forgot the imaginary unit i which can 
appear according to the sign of the term under the square root. 

In the large v limit, b\ is near 1 and b 2 is near y and the lowest relevant order 
for the potential for t around y is: 

-^v(y - 1) {{y - b 2 ) - ^y{y - l)x 2 ) 

which is the potential for a harmonic oscillator with energy (v/4h 2 )(y — l){y — b 2 ) 
and oscillator strength (v/4h 2 )y(y — l) 2 . This gives y — b 2 = 2h^J y / v{n 2 + 
1/2), which is exactly the semiclassical result at this order. Similarly, 1 — by — 
2H/^/v(nx + 1/2). In order to choose the terms which will contribute to the 
next order, we should remember that for the harmonic oscillator, x 2 is of order 
v ~ l l 2 . We therefore get corrections stemming from the variation of the oscillator 
strength, from the factor in front of y — b 2 and from anharmonic terms in x 4 
which stem in part from the one in r and in part from the product of the two x 2 
terms occurring in both factors t — b{. The order one corrected potential is: 

1 ;v ({y - b 2 ){y - 1 + -^L(m + ]-)) - \y{y - 1)V 



Ah 2 V Jv K 2 



2> 

With the anharmonic oscillator Schrodinger equation written as — \1>" + l/4(x 2 + 
Ax 4 )^ = Ety, first order perturbation theory yields the energy levels E n = n + 
1/2 + 3/4A(2n 2 + 2n + 1). This result is easily obtained writing x A in terms of 
creation and annihilation operators. Putting these corrections together, we get: 

, fy, 1, h 2 ( 2Jy , l w 1 1 1 „ \hy-\ 



v K 2' v\(y-ly 2 /v 2' 2 V 2' 8t/-l 

Note that this only differs from the semiclassical result by the last term, inde- 
pendent on the n.j. 

A similar computation yields the potential developed around t — 1 and the 
first order correction to the energy of the corresponding anharmonic oscillator 
gives 

1 - 61 = 2^-=(ni + -) + — ( 7 V lX (^i + ~){n 2 + -) - -(m + - N2 



V 2 y wV(y-l) v 2 yv 2 y 2 V 2 y 83/ -1 

Whatever happens outside of these two regions is exponentially suppressed for 
large v, since the wave function of the harmonic oscillator are already exponen- 
tially small there. Combining the values of b\ and b 2 , we obtain for the energy 
at this order in 1/v a result which differs from the semiclassical one of eq. Q15p 
only by a constant 



Eo — E\yxb ~ —ti 2 



This nicely reflects our findings on Figs. 2 and 3 where h = 1. 
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6 Conclusion. 



The interaction of the study of the Neumann model and the generalized Lame 
equation allowed us to give precise results on the behaviour of the spectrum in 
some limits. The generalized Lame equation has independent interest since it 
is the most general linear differential equation of second order with only regular 
singularities at finite distance. 

We were able to precisely count the different cases for the positions of the 
zeroes of the potential, showing in particular that they are all real. It is also 
remarkable that there is a constant error in the WKB approximation. This error 
is however independent on the level, so that it does not appear in transition 
energies. In this integrable case, the WKB approximation gives pretty good 
results even for low excitation numbers and the further study of its behaviour 
around turning points at low energy could make it even more useful. 
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